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ABSTRACT 

We have investigated the intrabinary shock emission from the redback mil¬ 
lisecond pulsar PSR J2129-0429 with XMM-Newton and Fermi. Orbital modula¬ 
tion in X-ray and UV can be clearly seen. Its X-ray modulation has a double-peak 
structure with a dip in between. The observed X-rays are non-thermal dominant 
which can be modeled by a power-law with P ~ 1.2. Intrabinary shock can be 
the origin of the observed X-rays. The UV light curve is resulted from the ellip¬ 
soidal modulation of the companion. Modeling the UV light curve prefers a large 
viewing angle. The heating effect of the UV light curve is found to be negligible 
which suggests the high energy radiation beam of PSR J2129-0429 does not di¬ 
rect toward its companion. On the other hand, no signihcant orbital modulation 
can be found in 7 —ray which suggests the majority of the 7 -rays come from the 
pulsar. 
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1. Introduction 

A new population of eclipsing millisecond pulsars (MSPs) has emerged in the last hve 
years, which is characterized with an orbital period Pb < 20 hrs and a companion mass 
Me ~ 0.1 — 0.5 Mq. They are commonly referred as “redbacks” (Roberts 2013; Hui 2014). 
Currently, 18 redback MSPs have been discovered. |§ Through coordinated multiwavelength 
searches, the population of redbacks is growing (Kong et al. 2012,2014; Hui et ah 2015). 
These MSPs play a crucial role in exploring the transition between rotation-powered system 
and accretion-powered system. 

To better understand such transition, it is important to probe the interactions between 
MSPs and their companions. In rotation-powered state, the collision between the pulsar 
wind and the mass outflow from the companion can produce intrabinary shock. The non- 
thermal X-rays from the accelerated particles in the shock region can be modulated at the 
orbital period. This has been observed in various redbacks/black-widows (e.g. Tam et al. 
2010; Huang et al. 2012; Hui et al. 2014). 

PSR J2129-0429 is one of the poorly studied redbacks, which was discovered in the radio 
pulsation search of the 7 —ray source 2FGL J2129.8-0428 with the Green Bank Telescope 
(Hessels et al. 2011; Ray et al. 2012). The dispersion measure suggests a distance of 
d ~ 0.9 kpc (Ray et al. 2012; Roberts 2014@). Its spin period and orbital period are 
Ps ~ 7.62 ms Pb ~ 0.64 d respectively (cf. Tab. 2 in Ray et al. 2012). Roberts (2014) have 
reported its surface dipolar magnetic held of R ~ 1.6 x 10® G. This implies a characteristic 
age and spin-down luminosity at the order of r ~ 4 x 10® yrs and R ~ 3 x lO®'^ erg s“® 
respectively. These suggest PSR J2129-0429 to be a young energetic MSP with a relatively 
high surface magnetic held (Roberts 2014). 

The minimum mass of its companion is found to be > 0.37 Mq by radio timing (Ray et al. 
2012 ). Optical observations suggest its companion is signihcantly bloated with a Roche lobe 
hlling factor of ~ 95% (Bellm et al. 2013). Its X-ray position of PSR J2129-0429 has recently 
been constrained by using the archival Swift XRT data (i.e. source J2129B in Linares 2014). 
Using the data obtained by XMM-Newton, Roberts (2014) has found signiheant X-ray orbital 
modulation from this system. However, the detailed emission properties of PSR J2129-0429 
have not yet been reported. In this Letter, we report the results from our investigation 
of PSR J2129-0429 by using the X-ray and UV data obtained by XMM-Newton. We also 
present the analysis of the 7 -ray data obtained by Fermi 7 -ray Space Telescope. 


®see http://apatruno.wordpress.com/about/millisecond-pulsar-catalogue/ for updated information. 
®http://fermi.gsfc.nasa.gov/science/mtgs/symposia/2014/program/10A_Roberts.pdf 
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2. XMM-Newton Observation &: Data Analysis 

The XMM-Newton observation of PSR J2129-0429 started on 28 October 2013 with a 
total exposnre of ~ 80 ks (ObsID: 0725070101; PI: Roberts). Using XMM Science Analysis 
Software (XMMSAS version 13.5.0), we rednced and hltered the data in a standard pro¬ 
cedure. The effective exposures for MOSl, MOS2 and PN after hltering are found to be 
~ 79.1 ks, ~ 79.0 ks and ~ 77.5 ks respectively. This X-ray observation covers ~ 1.4 orbital 
cycles of PSR J2129-0429. All the EPIC data are found to be unaffected by CCD pile-up. 

With all the EPIC data merged, we produced an X-ray image of the 3' x 3' held around 
PSR J2129-0429 ( see Figure H]). With the aid of the XMMSAS task edetect_chain, the 
X-ray position of PSR J2129-0429 determined by this observation is RA=21^29™45.250® 
Dec=-04°29'07.95” (J2000) with a statistical uncertainty of ~ O.l"0 This is consistent with 
the X-ray source J2129B reported by Linares (2014) but with the position better constrained. 

In Fig. [U we note that the X-ray image of PSR J2129-0429 is slightly extended toward 
west. Such feature can also be seen in the images from individual cameras. In MOSl/2 data, 
it is detected as a distinct source at RA=21^29™43.749® Dec=-04°29'07.83” (J2000) with a 
statistical error of ~ 0.7”. We do not End any identified object within 5” from this position 
in SIMBAD and NED. Also, we do not End any optical/IR counterpart in the USNO-Bl.O 
and 2MASS point source catalogs. Within 8” around its nominal position, we have collected 
~ 230 net counts from all cameras. We found that its spectra can be described by an absorbed 
power-law model with Ah < 5 x 10^° cm“^ and a photon index of P ~ 1.6. Its absorption- 
corrected Eux in 0.3 — 10 keV is fx ~ 1.5 x 10“^^ erg cm“^ s“^. Given the moderate angular 
resolution of XMM-Newton, we cannot unambiguously determine whether this feature is 
indeed a distinct source or a bow-shock nebula associated with PSR J2129-0429. Since the 
nature of this feature is uncertain, we excluded its contribution in all subsequent analyses. 

Before any temporal analysis performed, we applied barycentric correction to the arrival 
times of all the events by using the updated planetary ephemeris JPL DE405. We extracted 
the events from a circular source region of 15” radius centered at the X-ray position of 
PSR J2129-0429 so as to exclude the aforementioned feature in the west. After subtracting 
the background by sampling the events from the source-free regions in individual cameras, 
there are 718 cts, 732 cts and 2382 cts extracted from MOSl, MOS2 and PN CCDs, re¬ 
spectively. For improving the photon statistic, we merged all the EPIC data and folded the 
background-subtracted X-ray light curve at the orbital period of PSR J2129-0429, which is 
shown in the top panel of Figure [2l Since the radio timing model of PSR J2129-0429 is not 


Absolute astrometric accuracy of EPIC is ^ 1.2 . 
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publicly available, we arbitrarily chose the day that this observation started (i.e. MJD 56593) 
to be phase zero throughout our investigation. 

The maxima of the X-ray modulation are found at the orbital phases 0 ~ 0.1 and 
0 ~ 0.45. In between these two maxima, there is a dip at 0 ~ 0.25. We have inspected the 
background-subtracted light curves from individual cameras and found the dip at 0 ~ 0.25 in 
all three cameras. Therefore, we conhrmed that this feature is genuine. We further examined 
the X-ray orbital modulation of PSR J2129-0429 by dividing the EPIC data into soft band 
(0.3-2 keV) and hard band (2-10 keV) and investigated how does the X-ray hardness ratio, 
which is dehned as (hard-soft)/(hard-|-soft), vary across the orbit. We found that the X-ray 
hardness ratio decreases in 0 ~ 0.45 — 0.7 (see the middle panel in Fig. [2]) which is apparently 
in unison with the decrement of X-ray intensity (top panel of Fig. [2]). On the other hand, at 
phase 0 ~ 0.25 where a dip has been observed, there is an indication for the enhancement 
of X-ray hardness. 

The optical monitor (OM) on-board XMM-Newton has also observed PSR J2129-0429 
in fast mode with UVWl hlter (291 nm) for a total integration time of ~ 75 ks. PSR J2129- 
0429 is the only source detected in this OM observation. This provides us with temporal and 
photometric information of PSR J2129-0429 in UV regime. The background-subtracted UV 
light curve of PSR J2129-0429 was automatically extracted by the XMMSAS tool omfchain. 
After barycentric correction, we folded the UV light curve at the orbital period of PSR J2129- 
0429 with phase zero dehned at the same epoch for the X-ray light curve. The resultant UV 
light curve is shown in the bottom panel of Figure O Two UV peaks with a phase separation 
of ~ 0.5 are observed in one orbit. The minima at phases ~ 0.25 and ~ 0.75 apparently 
coincide with the dip and the minimum of the X-ray orbital modulation. For converting 
the magnitude into energy hux, we adopted a scale that the hux of Vega corresponds to 
0.025 mag in UVWl hlter. The energy hux changes from ~ 5.5 x 10“^^ erg cm“^ s“^ A~^ 
to ~ 1.0 X 10“^® erg cm“^ s“^ A~^ across the orbit. 

For investigating its X-ray spectral properties, we extracted the source and background 
spectra from the same regions adopted in the temporal analysis. The spectra obtained from 
all three cameras are htted simultaneously to the tested models. All the uncertainties quoted 
in this paper are la for 2 parameters of interest (i.e. Ax^= 2.3). 

We found that a simple absorbed power-law (PL) model can describe its phase-averaged 
X-ray spectra reasonably well {x^= 59.98 for 59 d.o.f). The best-ht model yields a column 
density of Ah < 6 x 10^® cm“^, a photon index of P = 1.25 ± 0.04 and a normalization of 
(1.84 ± 0.06) X 10“® photons keV“^ cm“^ s“^ at 1 keV. To investigate whether the phase- 
averaged spectrum requires an additional thermal component for modeling, we added a 
blackbody (BB) on top of the best-ht PL model. It results in a goodness-of-ht with 54.45 
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for 57 d.o.f. which indicates that the additional BB component is required at a conhdence 
level of ~ 94%. The PL+BB £t yields A^h < 3.9 x 10^° cm“^, P = l.ldlgQ®, a power-law 
model normalization of x 10“® photons keV“^ cm“^ s“^ at 1 keV, a blackbody 

temperature of kT = 0.16lg;Qg keV with an emission radius of i? = 92.9l4g^3^(iikpc m, where 
dikpc is the distance to the pulsar in units of 1 kpc. The best-fit PL-I-BB model and the 
observed spectra are shown in Figure [3l The unabsorbed energy flux in 0.3 — 10 keV is 
fx ^ 2.1 X 10“^^ erg cm“^ s“^. The BB component contributes ~ 2% of the total flux in 
this band. 

Motivated by the variation of X-ray hardness across the orbit (cf. Fig. [2]), we performed 
a phase-resolved analysis to investigate how does the emission nature vary with orbital phase. 
We divided the orbit into two intervals, 0 = 0.0 — 0.5 and (j) = 0.5 — 1, which encompass 
the peak and the trough of the orbital modulation respectively. For (j) = 0.0 — 0.5, we found 
that a single PL model is already sufficient for modeling the observed spectrum (x^= 57.57 
for 58 d.o.f) which yields Nu = 5. 2 ^ 5 * 2 ^ x 10^® cm“^, P = 1.10 ±0.06 and a power-law model 
normalization of 2 . 91 q 2 ^ 10“^ photons keV“^ cm“^ s“^ at 1 keV. Adding a BB component 
does not result in any improvement {x^= 57.51 for 56 d.o.f.). The best-fit PL±BB model 
in this phase interval yields A"h = 2 . 9 I 2 V ^ cm“^, P = 1.101 q 24, ^ power-law model 

normalization of 2.9lo;g ^ photons keV“^ cm“^ s“^ at 1 keV, kT = O.lblg;}® keV and 
R = 95.6]'^g5^/(iikpc m. Although the BB parameters are poorly constrained in 0 = 0.0 — 0.5, 
we note that their best-£t values are consistent with those inferred from the phase-averaged 
analysis. 

On the other hand, in 0 = 0.5 — 1, we found that the BB component is required at 
> 99% confidence level. In this phase interval, the PL±BB fit yields A^h < 4.5 x 10^° cm“^, 
P = 1.13lg26, ^ power-law model normalization of 7 .II 2 1 ^ 10“® photons keV“^ cm“^ s“^ 
at 1 keV, kT = keV and R = 84.4tltid ikpc m. The corresponding goodness-of- 

fit is 48.81 for 52 d.o.f. The BB parameters are fully consistent with those inferred 
in 0 = 0.0 — 0.5 and the phase-averaged analysis. This suggests the thermal component 
provides a constant contribution in all orbital phases. 

We have checked the robustness of the quoted spectral parameters by repeating all the 
aforementioned spectral fits with the background spectrum sampled from various source-free 
regions. Within la errors, the parameters inferred from independent fittings are consistent. 


6 


3. Fermi LAT gamma-ray observations 


Gamma-ray data were obtained, reduced and analyzed using the Fermi Science Tools 
package (v9r33p0), which is available from the Fermi Science Support Center 0. Events in 
the reprocessed Pass 7 “Source” class were selected and the P7REP_SOURCE_V15 version 
of the instrumental response functions were used. To reduce contamination from the Earth’s 
albedo, we excluded time intervals when the region-of-interest (ROI) was observed at zenith 
angles greater than 100° or when the rocking angle of the LAT was greater than 52°. 


For spectral analysis, we used photons between 0.1 and 300 GeV within a 21° x 21° ROI 
centered at the position of 2FGL J2129.8-0428 (Nolan et ah 2012). We performed binned 
likelihood analyzes with the gtlike tool. For source modeling, all 2FGL catalog sources (Nolan 
et ah 2012) within 12° of the ROI center, the galactic diffuse emission (glLiem_v05_revl.ht) 
and isotropic diffuse emission (iso_source_v05_revl.txt) were included. For sources more than 
10° away from the position of 2FGL J2129.8-0428, the spectral parameters were hxed to the 
catalog values. 


For spectral analysis, we used Fermi LAT data collected between 4 August 2008 and 17 
December 2014. We modeled 2FGL J2129.8-0428 with a simple power law 
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The ht with a simple power law gives P = 2.25 ± 0.06 and F^ = (1.57 ± 0.20) x 10“® 
photons cm“^ s“^ above 100 MeV, and a test-statistic (TS; Mattox et ah 1996) value 
238.2. The power-law spectrumis consistent with that reported in the 2FGL catalog (Nolan 
et ah, 2012). Although the reported 2FGL spectrum is not signihcantly curved, we also 
ht 2FGL J2129.8-0428 using a power law with exponential cutoff (PLE). This gives P = 
1.65 ± 0.21, Ec = 3.0 ± 1.1 GeV and E^ = (1.03 ± 0.21) x 10“® photons cm“^ s“^ above 
100 MeV, and a TS value 258.4. The likelihood ratio test gives 2A log(hkelihood) ~ 21.8; 
therefore the PLE model is preferred over the PL model at a statistical signihcance of ~ 4.6(T. 

We also constructed a long-term light curve with 3-month bins, and did not hud signif¬ 
icant step-like hux change like the one observed for PSR J1023-I-0038 in July 2012 when the 
latter changed from a radio MSP phase to a LMXB phase (e.g. Takata, et ah, 2014). 


^^http://fermi.gsfc.nasa.gov/ssc/data/analysis/software/ 
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We further examine for possible gamma-ray orbital modulation using 0.1 — 300 GeV 
photons obtained with Fermi observations over 6 years (4 August 2008 to 11 November 
2014). About 2860 counts were obtained in a circular region of 1° around 2FGL J2129.8- 
0428. With the same phase zero adopted in X-ray and optical investigations, the probability 
to randomly obtain a gamma-ray signal at least as significant as the one we obtained at 
the orbital period determined from optical observations, is 0.011 and 0.1 inspected through 
H statistics {H = 11.3) and test (y^=33.1 in 24 d.o.f.), respectively. So, we conclude 
no significant orbital modulation is detected in the gamma-ray band. The lack of orbital 
modulation does not support the shock origin of gamma-rays. A spectral cutoff at several 
GeV has been found from many gamma-ray MSPs (e.g., Abdo et al. 2009), and therefore 
the majority of the gamma-ray emission may come from the pulsar magnetosphere. 


4. Discussion 

We have investigated the X-ray and UV emission from PSR J2129-0429 with XMM- 
Newton. Both X-ray and UV light curves are signihcantly modulated at the orbital period. 
The X-ray emission is found to be non-thermal dominant with a small thermal contribution 
which presumably from the neutron star surface. While the intensity of the non-thermal 
emission varies across the orbit, we found the thermal component is essentially constant in 
all orbital phases. 

The X-ray modulation has a double-peak structure with a dip in between. This is similar 
to the original black widow PSR B1957-I-20 (Huang et al. 2012), in which the dip appears 
when the secondary is between the pulsar and the Earth and its X-ray modulation can be 
explained by the Doppler boosting of the shocked pulsar wind that wraps the secondary star. 
Since we do not have any information of the radial velocity of the companion of PSR J2129- 
0429 from either optical spectroscopy or radio timing observation, we could not determine the 
geometry of its orbit with the current data. However, the phase of the X-ray dip/minimum 
align with the phase of optical minima, suggesting the X-ray dip/minimum are at either 
inferior conjunction or superior conjunction. In this scenario, the X-ray modulation can be 
caused by the Doppler boosting effect (see Huang et al. 2012 for a detailed discussion). 

The optical modulation from the black widow and redback MSPs are usually observed 
with a single broad peak, since the hemisphere of the secondary star is heated by the irra¬ 
diation of the strong pulsar radiations. The double peaks with 0.5 phase separation seen 
in Figure [2] suggests that the optical modulation of PSR J2129-0429 is not caused by the 
irradiation of the strong pulsar radiation. Hence, it is expected that PSR J2129-0429 be¬ 
longs to a special case of the redback. Table 1 summaries the parameters of the redbacks 
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in the Galactic field (see Roberts 2012; Romani & Shaw 2011; Kong et ah 2012); the 
projected separation between two stars (asini, sixth column) is calculated from (Frank et 
ah 2002), asinf = 2.9 x ( Mjv 5 / Mq )^/^(1 + g)^/^(Po/day)^/^cm, where Mjvs = 1 . 4 Mq is 
the neutron star mass and q = Mmin/Mj^s with being the inferred minimum mass 

(fifth column) of the secondary and i is the orbital inclination. In all the following estima¬ 
tion/modeling, we assume a ~ asini. We further calculate the typical size of the Roche-lobe 
from r'L = 0.462a[g/(l + q)Y^^- The last column shows the ratio of the stellar luminosity 
without heating {Lstar) and the pulsar luminosity deposited onto the stellar surface {Lheat), 
that is ?7 = Lstar/Lheat- Siuce the lack of the optical observations for most of the redbacks and 
the secondary star is a non-degenerate star, we infer the stellar luminosity without heating 
using a simple relation Lgtar oc with (3 ~ 3.5 for a main sequence star. For example, the 
effective temperature of the unheated side of the secondary is Tef/fl ~ 2900K for PSR J2339- 
0533 (Romani & Shaw 2011) and ~ 5600K for J1023-I-0038 (Thorstensen & Armstrong 2005). 
The inferred stellar luminosities are Lgtar ~ ~ 1-8 x lO^^erg s“^ for J2339-0533 

and ~ 5.2 x lO^^erg s“^ for J1023-I-0038, respectively, which can be fitted by Lstar oc 
with a index (3 ~ 3.5. In the calculation, we used the parameter of PSR J2333-0533 for 
the normalization. The heating luminosity T/ieat is calculated from L/ieat = fnLgd, where 
= (1 - v/l-ri/a2)/2. In Table 1, we see that the stellar luminosity of the secondary 
of PSR J2129-0429 can significantly exceed the heating luminosity, suggesting the heating 
effect is negligible in the optical light curve. It is also interesting to note that the redback 
PSR J1723-2838 has similar orbital properties to PSR J2129-0429. PSR J1723-2838 may be 
another candidate that shows no heating effect in the optical light curve. Apart from the 
pulsar wind, the irradiation of the GeV gamma-rays from the pulsar is a possible heating 
mechanism of the secondary star in the black widows/redbacks (c.f. Takata et al. 2010, 
2012). In such a case, the absence of heating of PSR J2129-0429 suggest the gamma-ray 
beam does not direct toward its companion. 

Assuming that the companion of PSR J2129-0429 fills the Roche-lobe, we constrain the 
viewing angle of the system by modeling the UV light curve. We assume the Roche-lobe 
filling factor to be unity and the observed UV photons are produced on the Roche-lobe 
surface (i.e. optical depth is unity on the Roche-lobe surface). The Roche-equipotential is 
given by (Frank et al. 2002) 


^,,2 1 2q 1 X fX 2q \ 

^ix.y.z) = - , H- - - , -H-- — -h 

i+9y(f)" + (!P + (f)" i+9y(f-ip+(!P + (f)^ 9+iy 

(3) 

where we assume two stars orbit on x — y plane, and we used q ~ 0.2 of PSR J2129-0429. 
The Roche-lobe surface is defined by the plane ^{x,y,z) = <h(a;ii, 0, 0), where the inner 


Lagrangian point xli is calculated from 


a I'cs 
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1 1 

TT~qJ^ 


q _ 1 , ^ _ g 

1 + q — 1)^ a 1 + q 


0 , 0 < xli < a. 


( 4 ) 


We divided the Roche-lobe surface into many segments. For the limb darkening to each 
segment, we used a simple linear law, /(/r) oc 1 — /?(! — /i) (van Hamme 1993), where /i is 
the cosine of the angle between the local vector normal to the Roche-lobe surface and the 
Earth viewing angle. We chose a typical limb darkening coefficient k = 0.5. We ignored the 
effects of the gravity darkening (Lucy 1967) for simplicity. Hence the observed brightness 
of each segment is proportional to the intensity /(/r) times the projected segment area. In 
the bottom panel of Figure [21 we compare the UV light curve and the model light curves for 
various Earth viewing angles where ^ = 90° corresponds to the edge-on view. In Fig. [21 we 
can see that a larger viewing angle is preferred. This can explain the presence of extensive 
radio eclipse and the X-ray orbital modulation caused by the Doppler boosting effect. Since 
the current errors of the optical data are huge, high signal-to-noise optical observations will 
be required to tightly constrain the viewing geometry. 
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Fig. 1.— X-ray image of 3' x 3' field around PSR J2129-0429 as observed by XMM-Newton 
with all EPIC data merged. The black and white crosses illustrate the positions of the pulsar 
and a nearby feature respectively. 
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Fig. 2.— The background-subtracted light curve of PSR J2129-0429 as observed by XMM- 
Newton in 0.3 — 10 keV with the data from all EPIC cameras combined {upper panel), the 
X-ray hardness variation {middle panel) and the UV light curve as observed by XMM-Newton 
OM with the model light-curves overlaid {lower panel: see §4 for details). The orbital period 
adopted for folding is 0.64 day as reported by Ray et al. (2012). The epoch of phase zero is 
set at the MJD 56593. Two periods of orbital motion are shown for clarity. 
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Fig. 3.— The phase-averaged X-ray spectra of PSR J2129-0429 as observed by XMM-Newton 
PN (upper spectrum) and MOSl/2 cameras (lower spectra) and simultaneously htted to an 
absorbed power-law plus blackbody model (upper panel) and contribution to the y statistic 
(lower panel). 
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h 

J1023+0038 

1.7 

5 

0.2 

0.2 

1.2 

2.7 

0.8 

J1628-32 

3.2 

2 

0.21 

0.16 

1.2 

2.6 

1.2 

J1723-2837 

1.9 

5 

0.62 

0.4 

2.6 

7.2 

6.3 

J1816+4510 

3.2 

5 

0.36 

0.16 

1.7 

3.7 

0.4 

J2215+5135 

2.6 

6 

0.18 

0.22 

1.1 

2.5 

0.9 

J2129-0429 

7.6 

4 

0.63 

0.37 

2.6 

7.1 

6.4 

J2339-0533 

2.3 

2 

0.2 

0.075 

1.1 

1.9 

0.1 


Table 1: Parameters of redbacks. Pg] spin period. Lgd', spin down power. Po] orbital period. 
Mrain', inferred minimum mass, asini; the projected semi-major axis, size of the Roche- 
lobe. 7]] ratio of stellar luminosity without heating to heating luminosity. 



